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This work can Toe Toasically diTided in to two parts. 
In the first part, a method has Toeen developed for the 
evaluation of the Integral of Square of the Derivative 
of Error (ISDB) of time invariant linear system excited 
by deterministic input. The method is relatively simple* 
in theory and easi^w to compute. 

And in the second part the design of linear time 
invariant system has been formulated as a non-linear 
programming problem and solved by computer using well- 
knovm. gradient method, A weighted sum of different order 
of ISDE has been considered as a criterion to be 
optimised. 
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GHiPTER - I 


BYiiU.A!ri01T OP PERPORJIAFCE IMEGR^ 

1' IIT TRODUOI IQR; 

Oonsiflerable amount of ■'.'rork has been done on the 

eTaluetion of di-fferent types of integrals involving the 

1 2 

error of the linear time invariant system. Most of them ’ 
are ba^sed on time domain analysis and involve matrix transfer— 
m< 3 tions and Lyapunov’s functions. 

In this chapter a method has been developed for the 
eva.lur.tion of the integral of the following type 



which is the integral of the square of k th derivative of 
error, represented an ISDE(k). 

The method is based on Parseval’s Theorem and hence 
is a frequency doma,in analysis. The order of the derivative 
can be any thing even more than that of the system. One 
advantage of this m.ethod is thah once the value of ISDE(k) is 
evaluated the values of ISDE(i) for i = 0, 1, 2, . . , , 3&-1 
can easily be computed. 
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-p.^iSZVML'S 


Let Y(s) be the Laplace transform of 7(t) defined by 


Y(s) 


, 

/ y(t) e”^ dt Re [si > 


a 


.iiid let y(t) and Y(s) satisfy the following conditions 


lim 

± oa 


y (■fc) = 0 


and 


lijn 


S -♦ OO 


Y(s) = 0 


Then the follovring time and frequency domain relation of 


Parseval's theorem will be true 

CO 0“ 

/ y^(t) dt = f Y(s) Y(-s) ds (1.1) 

0 -jco 


where the region of convergence of the integrand in the right 
hand side of (1,1) centres along the imaginary axis in complex 
plane and is of width less than 2a. And hence the contour 
of integration for (1,1) can be taken along the imaginary axis 
and closing the entire left half plane. 

5* EYAL ITATION of ISLE : 


let X (s) be the Laplace transform of the error of the 


linear system given by 

E(s) ' 


X(x) = 


Q(s) 


m . 

r = o r 
n i 

i = 0 


m < n - 1 


(1.2) 
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Let Q(s) be the Hurwltz polynomial, 
loped for the evaluation of 


rSDB(k) 


I ii? 



dt 


The method will he deve 


... (1.5) 


where ISDE(lc) corresponds to the k th derivative of the 
error. 


5 . 1 Extensio ns of Parseval * s Theorem : 
Let y(t) 


d^ 


x(t) 


then 


■Y(s) 


dt 

L [y(t)] = L 


,k 




. k-i (l{>.i-.i) 

= s^X(s) - S 2:(o) 

i=o 


1 . 


mhere s is f'-e i th pov/er of s and 


(k-i-i ) 

X (o) = 


lim 


,k-i-.1 


t-O L ^^k-1-1 


X (t) 


Substituting the above, (1.1) gives 




. k-1 (k-i-l) 

s^X(s) - E s^ x(o) 
i=o 


. k-1 


(k-i-1) 


(-s)^ X(-s) - 2 (-s)^ X (o) I ds 

i=o • 
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jco 


Ic 


21c 


2 TCj - 3 = 


f . (_1) B X(s) X(-s) ds 


jco -j 


2%:. 


f 


^ (k-i-1 ) 

[ 2 (^1)^ X (o) .X(s) ] ds 


- 3 ^ x=o 


f I 2 (-1)^ s X (o) X (-s)|ds 

2113 - 3 - 


+ 


3°° 

r 


ic-i ]c-1 


;',T 3 


S 

3-0 


^ (k-.1-3) . k (k-1-i) ^ 

S s X (o) (-1) X (o) j ds 


1=0 


‘Bie oasni. integral, in the above expression wil.l 
vanish, because the integrand doesiiot have any pole inside 
the l.eft half pl^ane and hence it is analytic there. The 
fourth integral wi’’.! also vanish because of the following 
property : 


.1 

2713 


f s ds . = 1 if i = -1 and 

= 0 if 4 ¥ -1 

And hence ISDS(k) can be written as 


ISDE(k)r--— (-1)^ V X(s) X(~s) ds 

2 Ttj - 3 «= 


1 __ 

27i;3 


200 

f 


i ]^+-| (k-i-l) , 

s (- 1 ) 3 ^ X (o) X (s) j ds 

i=o 


V«- ( 1 . 4 ) 
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3.2 Spect_r^ PaGtorization^ ’ ^ 

Le-t PCs) be a rational function of complex variable 
s vtbicli b_as poles on both left and right half of the complex 
plane. Then P(s) can be seperated into two additive parts, 
one corresponding to the poles in the left half plane and 
the other to the right half plane. This process of separation 
is known as spectral factorization. 

Since in Parseval’s Theorem the integrand, which is 
of the form 

E(s) ~ X(s) X(-s) 

has poles in both half of the complex plane to evaluate the 
integral, the corresponding spectral factored term of the 
integrand will be needed. 

Tfith X(s) as given in (1,2), the spectral factored 
terms can be written as given below: 


X(s) X(-4 = “Ilf} 


P(-s) 


u(s) 

WJl 


+ 


IJlzSl 


( 1 . 5 ) 


where U(s) 


n.-1 

2 

i=o 


UpS 


i 


the order of P(s) is alv/ays (n-1 ) irrespective of the 
order of P(s). 

In the above expression U(s)/q(s) can be considered 
to be the sum of all partial fractions of X(s) X(-s) having 
poles on the left half plane. 
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'Ilie expression, in (1.5) can also be vnritten an 
P(s) P(-s) = U(s) Q(-s) + U(-s) Q(s) 


On equating the like powers of s from both sides of the 
above eouation a matrix equation of the form 

Rp = 23u ..,.(1.6) 

can easily be ob-i'a-ined. 

where R and S are ’ n x n matrices and 

jD, u a,re n x 1 coloumn vector. 

The vector u. will give the coefficients of UCs) in (1.5).- 
5 * 5 I5^.3i aples_: 

The follov/ing tv/o exajaples will illust.rate the 
techniaue involved in spectral factorisation. 


Example -1 


let 

X(s) = 


^7 


p. 


s+qgs' 


and 


U(s,) _ 

7(77 ■ 


u^+ u. s 

o 1 


q.o+q^s + q2 


Hence 

X(s)X(-s) 




+ qgS 


u^ + u^s 


Uo - s 




4 * 
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G ro s s-mtil t ipl ic at io n gives 

2 

Pq'’ ^^0 **■ s + q2 s ) 

+ (u^~u^s) (q^+q-i® + <l2 


Oollecting the ooeffioients of like powers of s following 
two 8qun‘;:ions are obtreined; 



0 

s : 

!l 

CM 

O 





2' 

s : 

0 

2l2 '^0 “ 

1 


(Thus 

i:he matrix 

eqr'at i'-n; 

L given in (1 

.6) can 

be written as 


[Po °i 

; i 

Po'l 

2 fa 0 1 

1 0 

1 ■ I 

IN 

i 

[ 


[o 0 J 

1 0 J " 

Ua - 

Ih 

I 

ibcarij 

Let 

ole 2 

.X(s) = 

P(s) 

_ Po + PiS ■ 

, 2 
■b Pg S 


T(s) 

2 

Iq, + q^s -Hq^s + q, 

“3 

5® 


In the similar way as mentioned in example 1', expression (1.6) 
can be derived. In this case 


Po 

• 1 

0 

0 ■■ 

O 

o 

O 


^ol 


'u 

0 

£ = jPq 

i 


^0 

1^12 ^0 1 

E = 

1 

f i 

: 1 

u = 

Pi 

|o 

0 

Pg: 

1 i ’ : 

1 I 0 ^2 

J 1 

i 1 

fV) 

1 


P2 
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Cm) 

5*4- E~yai -ua-t ion _of x(o ) 

■Tlie principal part of the Laurent Series expansion 
of X(s) is given by . 


X(s) = 


P(s) 

Q(s) 


where d = — — 
^ 2txd 


n=1 


r x(s) 

Jo as (1.7) 

S 


This expansion is va.lid out side the circle centered 
at the origin, ifiaose radius is given by 

r. > max j z. | 

^ ■ i ^ 

where is the vector in the complex plane representing 
the ith. zero of Q(s), The contour of integration c can 
be ta,ken any where out side this circle. 

Es'pansion (1,7) can easily be obtained by long hand 
division of the numerator P(s) hy Q(s), Taking the inverse 
Laplace transform, (1.7) gives 

... ( 1 . 8 ) 


(1.9) 


x(t) = 2 

n=(5 




n 


n 


Prom (1.8) it is easy to see that 


(m) , 

^ (o) » <3^+1 


lim r d 


m 


laiJ" 


X(t)! 



9 


5»5 B\'‘aluatloH "Of ■ X(s) ds 

The above ban be very easily obtained by expanding 
X'(s) in terms of the Iiaureht Series given by (1.7). Integrat- 
ing the expansion term by term, and noting that 

s^ ds = 0 if ± ^ -.1 and , 

= 1 if i = . 

The following can be obtained 

2^5*^ '^Q ^ ' X ( s ) ds ^ ^^ 4*1 ^ io ) 

5.6 Evaluation of s^X(b) X(-s) ds 

The a,bove integral can be represented by 

m 4 ^ 4 S(f} 

- 4i4 ^ 

The second integral will vanish, if the bontour of integration 
is taloen in the left half plane, where the integrand is 
analytic. And hence 



Lei: the Iiaurent Series expansion around, the origin of IJ(s)/q(s) 
be given bj 


FM = I n 

Q(s) n*1 gU 



1 r u(s) 
^ -^c q(s7 


gn-l 


ds 


( 1 . 11 ) 


And hence by the sarae way as mentioned in section 5.5 


JL 

27tj 


r X(s) X(-s), ds = C, 


'm+l 


... ( 1 . 12 ) 


5,7 Value o;f ISDE(k) 

# 

The expressions so far "derived can now be used to. 
evaluate the expression for ISDE(k) given in (1,4). Thus 


ISDS(o) 


ISDE(k) 


~ C. 


= (-lVc 


2k+1 


■ k-1 
i*o 


^k-l ‘^k+i+1 


... ( 1 . 15 ) 

for k>_ 1 


where O's and d's are given by (1.11) and (1.7). 
Thus, 

ISr)l](o) « : 


ISfSd) s “ <^1 *^2 

ISi)E(2) = + d., - dg d^ 

iSBSCf ) * -0.^ ’ - d.^ , dg' + d^ d^ - ,..d^ ■ 



GHiPTER - II 


fopMjAIIoi of the PROBIM 


1. ISU^ROjUGTIOB 


Tlie classical method of S37'stem design is based on 
trial and error technique and depends very much on the 
er.perienc8 of the designer. , The possibi“!.it 3 r of formulating 
the design procodure o.s a non-linear programming problem has 
been mentioned b^y ITex^'toii^, But the appearance of the actual 


formulation in the liter.ature is relatively recent 


. 5,6 


The fortiiulation presented here is a more general one. 
In t ‘is rrethod the error of the system is expressed in terms 
of the laplace transforir}, foimij, as a ratio of two polynonials, 
and the weighted sum of the ISDE(lv) for different values 
of k, is minirai.?.ed by selecting proper valu.es of the co-effi- 
cients of the above two polynomials. The o.rder o.f the 
denominator polynomial should chirr be given and the numerator 
will automatically be fixed up by the me'^rhod. 
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2 * PERFORMANCE IM)EX ; 

For the opthTi-um design of linear systems many types 

of performance indices have been developed, Unfoirbunately 
are 

none of them ~i— i suitable and very general for the design, 

A detailed discussion of the different types of performance 

7 8 

indices ha.ve been discussed in many literature ^ . 

In this design method the system error is represented 

in its Iiaplace transform, form as 

m i 

• ^ ^i ® 

X(s) = E(s) = — m <_ n - 1 

n 

S a . s^ 
i=o ^ 

with e(t) an its time function. The following perlTormance 
index is proposed here 


1 = 


M 

2 

k=.o 


w. 


k 


ISPE(lc) 


M 

= 2 
k=o 




.k 


01 dt 


k 



( 2 . 1 ) 


The above can be found out by the method presented in 
Chapter 1. 

Suitable vreights can be selected to give proper 
stress to the desired derivative of the grror. The vadue 
of M depends on the choice of the designer. 
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3. GONSTBAIUTS; 

i) IThe derivation of the performance index y given in 

Chapter 1, shows that (2.1) is a function of C*s and d's 
which in turn a.re functions of £, £ and u, where the 
coloumn vectors H constitute the co-efficients 

of the numerator and denominator of the haplace transform 
of the error and thaf of the numerator of the spectral 
factored term, respectively, and is given by (1.2) and (1.5). 

Thus performance- index (2.1)can be written as 

J = f (p_,- £, u) .... ... (2.2) 

where £ is an n- vector, £ is an (n+1) vector and u an 
n-vector where as f is a scalar function, (n being the 
order of the system). 

In order to define u in (2.2) the expression given in 

(1.6) will have to be satisfied. To take care of this, 

(1.6) will be considered as a set of equality constraint 
on the design method, and v/ill be represented in vector 
matrix form 

^ 3s,> u) ~ 5. *** *•* ^^•5) 

where the function h is an n-vector. 

ii) The basic requirement of the system is its stability. 

Tha epjyQU of -fihe eystili is given by 
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m 

E(s) = — m < n - 1 

21 ^ 

2 Ci.s^ 

i»o ^ 

To see that E(s) he stable and minim'um phs-se type both 
the denomimtor and iiumerator polynomial must satisfy the 
South criteria for stability. 

To consider this, the first column of the Routh 
table will be evaluated, corresponding to the numerator 
and denominehor polynomials and set them to be grea,ter 
than zero. These constraints can be defined as 


£l > ° 

r^ = n vector 

E2 (e) > 2 

£2 == (n+1) vector 


iii) In the minimization procedure, to avoid the 
trivial solution, p will be set equal to unity and 

m 

S ^ 1 ••• ••• (2.5) 

iso *” 

The inequality constraints given by (2.4) and (2,5) can 
all be accomodated and represented in the vector matrix 
form 

£ (e> a) i 2. ••• ••• (2-6 r) 
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wl'-iela rn.a.v include aj^y other extra Inequality constraints 
req;aired for the design, and accordingly the order, of the 
oolo-umn vector g will change* 


4. IFE PROBL Ua 

Ihe design of tie linear system can now be mentioned 
.aS belo.w; 


rind out the suitable values for £'s' and ^*s 
of "bhe Iia,place transform of the error given by 


E(s) 


s 


z 

i =:0 



m < nr-1 


so that the performance index 


■ J ^ (e* 3.9 

is minimized subje.ct to the constraints 

B ( b » 2 .) 2 : 2 
^ h), * 2 

p <_ n vector, u n vector, £ = (n+1) vector, 
h = n vector-. 


where 



* 

Kie aboye is bhs general form of the ITon-llnear 
programming probleip and can be solved bj?' any of the well 
IcnoTOi teohn.icaieai Here in thid thesis the solution has 
been carried out by the Gradient Method, 

1H 5 SO LITT'IOCT; 

To take care of the ineqUo3.ity constraints, these 
W3111 be added to the cost function in such a wa,y that if any 
solution point does not violale a constraint, then tha,t cort- 
straint will be set equal to zero and hence its contribution 
to the cost function v/ill be nil, on the other hand if it 
violates any constraint then the square of the violated 
constraint vll]. be aided to the cost function with a high 
penalty attached to it. Thus the new cost function can 
be represented as 

f (e» It a) + (e» a) ••• (2.7) 

Ic X 

where j= 0 if (p, £_) 2 and 

= 1 if (E, e) < 0 

and = a small positive quantity. 

For a particular value of the minimization 
procedure is canried out, and stopped if the final solution 
does not violate any of the inequality constraints. But 
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if the solution violates any constraint then the value of 

is reduced and the computation is repeated. It has been 

9 

sho\TO that the proper solution will be obirained for a 
sequence of values of By. tending to zero-.- 

So far as the equality constraints are considered they 
are al v^ays kept satisfied during computation. To start 
the computation, some arbitrary values of £ are assumed 

and tbs constraints h (£, £, u) = 0 are solved for u. 

Then £, £ are modified by the gradient method and modified 
u is obtained by solving the equality constraint, and the 
process is repeated, 

S'* SEEilCTiniT 01' STEP SIZE I N TH E GRfflIBM'. METHOD 

To find out the proper step size in gradient method 
is a difficult problem. An improper step size may lead to 
instability or a large number of iterations and hence lo^g 
computation time, 

A method for selecting the step size is proposed here. 
This iTBthod will lead to the optimum point from the assumed 
initial point in the least number of iterations. 

In gradient method, at any initial point the 
gradient of the objective function is found out, and moved 
along that gradient upto a distance equal to the step size, 
and then gradient is found out at that point and the process 
is repeated untill the optimum is reached. 
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Figure 2,1 shov/s the contour of the function to be 
ninimized, in the ’co-ordinate space'. The gradient of the 

function at a^ny point 
is normal to the contour 
of the fu.nction through 
that point. This gradient 
is also tangent to some 
contour of lower walue. 

The value of the contour 
to which the gradient is 
tangent is the least 
value of the function that can be obtained by moving along 
this gradient. The distance of the point of tangency from 
the initial point is the optimum distance along this 
gradient which v,dll give the manimum reduction in the cost 
function. This length will be the step size for the 
current iteration. 

So the step sise at every iteration will be found 
out by moving along the gradient as long as the value of 
the cost function decreases and the constraints are not 
violated. The point at which the cost function will 
start increasing o^ the solution will violate any constraints 
will be the initial point for the next iteration i.e, a ne\¥ 
gradient will be found out at this point and the process 
will be repeated. 




CHAPTER ■- Ill , 

EXAI'IPI'ES ■AM RESUI.rS 


• ■ S Tai UlTIOH np ISD E(k) 

Let x(s) = 

s^5so+.6 

Eor this X(s) the equation given, in (1.6) 'becomes 


t ol t 

2 

’ 6^. ,0:. 


O 

■ 1 1 


'i 



iO -21 12 

S 1 

t 1 

_ 1 -5 : 




Solving the above, u is obtained as- 
ju‘ [2.08 V 

j =1 

I uj- [0.816 _ 

and hence spectral faotored term corresponding to the 
poles in the le-ft half plane is given by- 



0 .816 s + ,2 ..0 8 
s^ + 5 s + 6 


Thus 


= 

0.8167 

“4 ' 

- 13.50 

II 

0 

296.1 

02 = 

- 2.00 

O 5 = 

36,90 

0 

II 

5.1 

% = 

-103.50 


2.00 

dg * 

- 5.00 

d^ - 

13.00 

"4 = 

-35.00 

d5 = 

97.00 


= -275.00 
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a.nd hence 

ISDS(o) = = 0.8167 

ISDS(1) = -,0^ _ d.j = 4.90 

ISDE(2) = d. ~ d d = 3 9 

ISDE(5)' = -G^ - d^ dg + d2 d^ - d^ « 224.0 

[Dne above result can easily be very.fied by the following 
direct conputation 

X(s) = 2ft.V5 ^ _1 1 

s + 5s + 6 s + 2 s + 3 


Therefore x(t) = 


+ e 


•3t 


and hence 


n(t) = ^2e"^'’^ _ 3e~'^'^ 


x(t) = 4e“^'^ + 9e''^’^ 


x(t) = -8e 


-2t 


- 27e' 


-3t 


} f + be’“2^ ’ 


dt = 


b^ 

^ " 2^2 


+ 


2ab 

a^4- 


Therefore ISDE(o) — y x^dt — ^ ^ ^ = 0.817 

0 

eo 

ISDE(1) = /^ x dt = ^ = 4.9 

ISDE(2) = 7 x^dt = i|. + ^ + 72 ^ 
ISDE(3) ^ / X dt = gI-27 + 16.27 , 

o 465 
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2 • DBSI&I1 0? A UMT MJICBBATOR SYSTWi 


This article presents the design of a unit nunierato 2 
third order system. The bloch diagram configuration is 
■ showm in Tigure 5.1. 



Figure 5.1 


let the system transfer function be given by 


&(s) 


1 



s-^ + a^s + a2 s + 1 


suitable values of a^ and a 2 are needed so that the 
following is minimized: 


J 


M 

E 

1g=o 


V;, / 

k -b 


” ^d^’ 


L dt 


k 


e(t) 


7 2 


dt 


... (3.1) 


where e(t) is the inverse Laplace troaisforin of E(s), 
For the configuration shown in Figure 5.1, 


E(s) 

E(s) 

R(s7 


R(s) ~ C(s) 


0(s) 

R(s) 


= 1 - G-(s) 


s + a^s + ag s 


s + a^ s 4 - a 2 S + 1 
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Hence ;f03? a, step input of unit amplitude 

2 

. . s + a^s + 

E(s) = 

3 2 

s + a^ s + agS + 1 

Comparing the co-efficients with 

m 

S p^s^ 

E(s) « m < n - 1 ... (3.2) 

j Q.gi 
1=0 

following additional constraints are obtained 
% = 1 Po = ^1 P2 = f 13 = '• P-t = ^2 

The results of the optimum system for different values 
of M and = 1 for all i, are presented in Table 3.1. 

TABLE 3.1 


■ M 

J 

ai 

^2 

0 

1.5366 

1.3450 

1.7995 

2 

2.1576 

2.2589 

2.1608 

4 

4.1564 

2.0670 

2.0921 


It is to be noted here that the values of and 

7 

a.2 obtained for M - 0 is the same as reported earlier. 
The transient response of the optimum system for an unit 
step input is given in Eig. 3.3. 
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5* PSSIG-F Qj? A QEIIERAL THIRD ORDER SIS TM 


Dor the similar configuration as given in Dig. 3.1 
a general third order system can be designed. 


Let G(s) 



Let - bg , then for a unit step input 
^ s^ + (a^ - b^) s + (a 2 - ) 

s^+a^s +a2® 

Gomparing this with (3.2) the additional constraints 
obtained are 

= 1 and P 2 = 1 

The values of a*s and b's axe given by 

3 -} ^ 0-2 ^2 ~ *^1 ^3 ~ ^0 

= Ig- P, I’l = <11 - Pa 1=2 = % 

The optimum values of the a's and b’s are given 
in Table 3.2 for different M and with = 1 for all i. 

The normalized transient response for a step 


input is given in Dig. 3.4- 
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TABLE ~ 3.2 



0 0,116 4.30 4.25 2.05 4.30 4.-T0 2.05 

1 1,000 1.00 2.59 1.65 3.60 4.25 i;65 

2 1.155 0.58 2.05 1.65 5.50 3.96 1*65 


the BESIGCT 0? CO]\(gEESATOR 

Tlie design •becHnique developed in chapter two caii also 
be extended for tbe design of a compensator of a control 
system, TLe block diagram .of the s 3 rstem is given in Eig.3.2 
where is the given plant transfer function, and Oq(s) 

is the compensator transfer function which is required to be 
designed. 



FIGURE 3.2 


let Gp(s) = 


m 

S 

i=o 



n 

S 

i=o 


a. s 
r 


i 


m <_ n and = 0 


(3.3) 
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I'Toto ■nlia.t = 0 will, le needed here for the design, this 
condition is generally satisfied in most of the phnrsical 
•Diarrfea, 


'}i 6 r 0 f oro 


&(s) - G^(s) Q-^(s) = 


v 

S d 

1 . 

T 

n-S" 

/ 

k 


2 

G. 3^ 

i==o 

2 _ 


r 


E u-s 

S 

1 

o 

II 

Hi 

i=o 


< k 


... (3.4) 


iKs) _ 
R(sT 


1 + G ( 3 ) 


i 

E 0 , . s 
i=o 

k 

E 

i^o 


n 

y 

ft V. , ' 

Zj 

i-0 

? -i 

S , a. a 

k 

2 

•1 

O 

11 

i=o 


k , . 

S C.'S 
i=o ^ 


ni 


r=o 


X / - i 

% d|s^ 
iso ^ 




<5j5) 


ind hence .for a unit step inp'at 

n+k-1 - . 


E(s) 


S p.s^ 

i^O 


n+K 


(3.6) 


S q_|S- 
i=o " 


Thus in this case the a,dditional const rairct is that 

n+k~1 


n+k 

S 

. 1=0 


S l.^s 


1 


S p.s- 

i=0 


should have 


m . 
i=o ^ 
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as ' factor, ijoi remainder -after dividing hj S 

slioriid be ern„ia3- to zero or lose than sorae small preassigned 
tolerance value. Any constraint equivalent to this cs^n 
also be usod^ 


It is aJ-SO desirable that the innut excitation 
does T:ot saturate the plaiot'. hor thisj follov?ing constraint 
is posed on i(t) (sliouii. in Fig, 3.2), the inverse Lap^ce 
i: r ans f o rm of I ( s ) ': 



oo 



idt) 


clt 



(5^7) 


Ilie above can be evaluated by the help of Parseval*s 
Iheorsm l-mov/ing the transform of iCt). 
lou I(s) ^ E(s) l^(s) 


n+k- 1 
S 

t=o 


p.s 


1=0 


a,s^ 


n+lc 

S 

i-o 


k 




4. a 




C.s- 


To evaluate the left hand side of (3,7) G-^(s) will have to 

be found out 'with the help of (3.5) and (3.6), 

The following section gives an example of the 
compensator design problem. 



5-. AM'- EX.aiELE OE COMEElTS.iTOR DESIGN: 


Let 


,(s) = 


4 


1 


8 .]ld 


G- (s) 
0 


s(s/2 + 1) 
d^B 4- cIq 


0.1,25 s + |25s 


n o 

_■» ' • V</ 


■Tlie optiimam vad-ues of d's and C's are needed so that 
(■’5,0) Is minimised satisfying (3.-7). 


G-(s ) — 


d^Cs) CxpCs) 


d^e + 


■125 a^s"’ t (.125 C^t .250,^)s + .25C^ 


i i1250,j.s''.±, (v125Cq + ;-25Cij)s^+ 4?5 C^s__ 


E(£) 

.TtCs) 


1 + &(s) . 1250, r0+( ,- 1 25C^^+.:250,j )s^+(v25CQ-i-d:5 )s4-d 


fon.s for 'unit step 


-^(s) 




fi25G,3^ + (.1250^, + ,25C^),s + .250^ 


1250^3^ +(_,125C^+^.-250.>s^ + i,25Q^ + d|) s + d^ 


Thus the a,dditional Qoiistra.ints cs.n be vn:itten as 


T). 


^3 


= *-1250,^ + .250,, o; 


'B T) 

' 0 


= 2 qp - 4 


.■aid hence -the parameters of the oompensator is given b 3 !' 
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T}ae optim-uun values of c's aiid d's for M = 1 and , 2 and 
for Vi7j^t= 1 .for all i is given in fable 


1AJ3L1 3.3 


M J 

„ 1 

rl 

0 





0,14 

59.56 

27^75 

V “' - ■ 

' 10.76 

8 

2 5)41 

0^12 

59.49 

-2,02 

150^04 

8 

lie no realised tra: 

given in ?igi 3 , 5 . 

nsienb : 

rosponae- 

for a unit 

step 

is 


i 

J 


1.4 


1.2 


1.0 


,8 




A \ 

/ \ 






.6 


.4 


.2 


' r\ 


M=1 


I ; 




d.s + 

1 0 


GiS + 


,125s +.^5s 


! // 


PI&UKE 5.5 


Seconds 



GOEgLug^oiis 


A Giomple‘i:e nusierical method has been presented .'for 
the optimijin design of linea.r time invarifant systeml The 
acLyphtagb of the niettod is i-hat it doss not need siij 
exper^^ce of; the designer. Iny arbitrajojr choice of initial 
pa7^^]leter of the system can be made to converge td; the 
finaJ- optimal vaJ-'oe. 


A xien performanoe measure has been pro.posed for 
the optimum, design of linear s;/stein. .Prom the transient 
responhoe given in Fig, o 3? 3.4- end 3.5 it is obsc-rvod that 


the orfterion considered is a reasonable one^. YiTith proper 
cfcoic'e of'' -meights for the ISDECIc) transient response can 
be given p:cop6-r shape,; 


.iilthough the design procedure consists of large 
number i^f comput a.tions at ever-,?- iterations^ this is not 
■a problem no'w a. days because of the availability of high 
speed domputers. 


It has not been possible to insett a suitable , 
si.nj.ple criterion which will assure that the optiinum system 
■will be Qont roll able* 

As a fu-rther research in this subject following two 
important problems are prese.nted. 

1 ) liven two pol 3 /nomicals and E find out 

a simiple relation involving a's and b’s so that the 
above polynomial. s will not li.ave any ooimnon root. 



2 ) 


value 

used 


A detailed study 
s, botli negative and 
ill the cost function 


01 L-lie cost function fon different 
positive, of w^, the weights 
• in so find out a method to 


.^Lcct optimuiii vdJ-ues 


w. 

1 


pojrticular value of M, 
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